Dynamically Induced Locking and Unlocking Transitions in Driven Layered Systems 
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Using numerical simulations, we examine a simple model of two or more coupled one-dimensional 
channels of driven particles with repulsive interactions in the presence of quenched disorder. We find 
that this model exhibits a remarkably rich variety of dynamical behavior as a function of the strength 
of the quenched disorder, coupling between channels, and external drive. For weaker disorder, 
the channels depin in a single step. For two channels we find dynamically induced decoupling 
transitions that result in coexisting pinned and moving phases as well as moving decoupled phases 
where particles in both channels move at different average velocities and slide past one another. 
Decoupling can also be induced by changing the relative strength of the disorder in neighboring 
channels. At higher drives, we observe a dynamical recoupling or locking transition into a state with 
no relative motion between the channels. This recoupling produces unusual velocity-force signatures, 
including negative differential conductivity. The depinning threshold shows distinct changes near the 
decoupling and coupling transitions and exhibits a peak effect phenomenon of the type that has been 
associated with transitions from elastic to plastic fiow in other systems. We map several dynamic 
phase diagrams showing the coupling-decoupling transitions and the regions in which hysteresis 
occurs. We also examine the coexistence regime for channels with unequal amounts of quenched 
disorder. For multiple channels, multiple coupling and decoupling transitions can occur; however, 
many of the general features found for the two channel system are still present. Our results should be 
relevant to depinning in layered geometries in systems such as vortices in layered or nanostructured 
superconductors and Wigner or colloidal particles confined in nano-channels; they are also relevant 
to the general understating of plastic flow. 

PACS numbers: 83.50.-v,81.40.Lm,62.20.fq 
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I. INTRODUCTION 

A collection of interacting particles on an ordered or 
disordered substrate undergoes a depinning transition 
under an applied drive from a pinned state with immo- 
bile particles to an elastic or plastic sliding state |l|-[lOl|. 
In an elastic sliding state, each particle keeps the same 
neighbors over time [l-Q , while in a plastic slidingstate, 
the particles do not keep the same neighbors d, |3, Hol - 
Il4l |. Many systems show a transition from elastic depin- 
ning for weak substrates to plastic depinning for strong 
substrates that is associated with changes in the veloc- 
ity force curves and the depinning threshold as well as 
the onset of strong hysteretic effects [3, [III, ^^ [la] ■ For 
example, the peak effect phenomenon for the transport 
of superconducting vortices, where a peak in the critical 
depinning force or critical current appears as a function 
of temperature or magnetic field, has been associated 
with a transition from elastic to plastic depinning and 
is known to be related to the disordering of the vortex 
lattice [5|, [il| . Although the peak effect is primarily as- 
sociated with superconductors, it can in principle appear 
in other systems that undergo a transition from elastic to 
plastic depinning. Several recent studies of colloidal sys- 
tems and frictional systems have revealed peak effect-like 
behavior [i,[il[i^. 

Plastic depinning transitions exhibit a wide variety 
of characteristics depending on whether the substrate is 
random Q or periodic 0,11, [13, [l7l-[20|. The resuhing 
plastic flow ranges from motion in winding channels to 



avalanche behaviors with strong fluctuations, coexistence 
of large pinned and flowing regions, or transitions from 
two-dimensional (2D) mixing phases to one-dimensional 
(ID) decoupled channel phases. For higher drives, many 
systems undergo a transition from plastic flow to a more 
ordered flow state, such as a moving smectic state where 
the particles organize into ID channels that can slide 
past one another [l4 j_ la_ 2l| or partially ordered moving 
crystal states [] 
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Due to the complexity of plastic depinning phenom- 
ena in 2D and three-dimensional (3D) systems, numer- 
ous simpler systems have been proposed that still exhibit 
plastic depinning, such as particles flowing in a series 
of simple coupled layers. Many systems with depinning 
transitions can be modeled effectively as layered systems, 
including sliding charge density waves I ill 2 211 and vortices 
in strongly layered superconductors j23l |24| . Advances 
in fabrication techniques have made it possible to create 
nanostructured systems in which particles such as vor- 
tices move through effectively ID coupled channels [25| . 
Marchetti et al. considered a mean field anisotropic slip 
model of coupled channels of particles oriented parallel 
to a driving force [2g. Along each channel, the parti- 
cle interactions are elastic, but slip can occur between 
neighboring channels. In this model, when the depin- 
ning is elastic the channels show no slip or hysteresis, 
but for plastic depinning there is both slip and hysteresis. 
Further theoretical work on systems with only two lay- 
ers demonstrated that a transition from non-hysteretic 
elastic depinning to hysteretic plastic depinning occurs 



as a function of disorder [27|- For 3D layered sliding 
charge density wave systems, a transition from an elastic 
depinning to plastic hysteretic depinning has also been 
predicted [2^, along with a second coupling transition 
at higher drives when the charge density waves begin to 
flow coherently. Numerical work on 3D layered super- 
conducting vortex systems has shown a similar disorder- 
induced transition from elastic to plastic flow marked by 
a large increase in the depinning threshold when the lay- 
ers decouple, which is typical of a peak effect behavior 
[23, 24, 28]. We note that although the transition to plas- 
tic depinning is marked by an increase in the depinning 
threshold, in real superconductors the depinning thresh- 
old decreases as the superconducting critical temperature 
or field is approached due to the changing penetration 
depth. This feature is generally not included in most 
simulations. 

In this work we consider a simple depinning model con- 
sisting of two or more coupled ID channels of repulsively 
interacting particles, all of which are subjected to a driv- 
ing force and quenched disorder. This system could be 
experimentally realized using colloidal particles in cou- 
pled ID channels with a disordered substrate and driven 
by an external electric field [2y-\Bzl, coupled ID wires 
containing Wigner crystal phases 33 36], or vortices on 
corrugated substrates driven parallel to the corrugation. 
Our simulations of this model include much more detail 
than can be captured in mean field studies [2^ [23, Isjj ■ 
In addition to confirming many of the predictions from 
the theoretical studies, we find a rich variety of new fea- 
tures. We induce decoupling transitions by changing the 
relative strength of the disorder from channel to chan- 
nel, something that has not been considered in previous 
work on the depinning of layered systems. We also study 
commensuration effects by varying the particle density in 
one channel relative to the density in other channels. We 
drive all of the layers of particles; previous studies of this 
type of model considered transformer geometries with the 
drive applied to only one layer [3^, 123 • Despite the ap- 
parent simplicity of our model, we find that even the 
two layer system has a wide variety of dynamical phases 
and exhibits all the salient features found for elastic and 
plastic depinning phenomena, including a peak effect at 
the transition between the two types of depinning. We 
also show that the transport signatures of the dynami- 
cal coupling and recoupling transitions can be enhanced 
in systems where the channels do not all have the same 
amount of disorder, suggesting that experimental realiza- 
tions of this geometry may be an excellent way to probe 
these dynamical states. 

The paper is organized as follows. In Section II we 
describe our model and simulation technique. Section 
III discusses decoupling transitions in a two channel sys- 
tem as the spacing between channels is varied, as well 
as a peak effect appearing at the transition from elas- 
tic to plastic depinning. In Section IV, we consider two 
channel samples in which the strength of the disorder dif- 
fers in the two channels, where we observe drag-induced 
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FIG. 1: Schematic of the system. The particles (filled circles) 
are constrained to move along ID channels that are separated 
by a distance d. Each channel contains Nj particles as well as 
rij randomly placed pinning sites (open circles) that each have 
a maximum pinning force of F?. A uniform drive F^ is ap- 
plied to all particles in the positive x direction. The particles 
interact via a repulsive Yukawa potential with other particles 
in the same channels and in neighboring channels, (a) A two 
channel system with A^i = A'^2 and ni — n2. (b) Two channels 
with an equal number of particles in each channel, Ni — N2, 
but with more pinning in the upper channel, ni > 722. (c) 
Two channels with equal pinning in each channel, ni — n2, 
but more particles in the upper channel, A'^i > N2. (d) Mul- 
tiple channels each containing the same number of particles, 



iV, 



A^i, and pins, Uj = ni. 



pinning as well as hysteretic velocity-force signatures. 
In Section V, we find commensuration effects when the 
number of particles in each channel is varied for a two 
channel system, and we also study the effect of changing 
the particle-particle interaction prefactor. We introduce 
an eight channel system in Section VI and show that as 
disorder strength and channel coupling is varied, the be- 
havior resembles that of the two channel system. We also 
study representative examples of eight channel systems 
in which not all channels contain the same number of 
particles. We conclude with a summary in Section VII. 



II. SIMULATION 

We model an array of M coupled ID channels sepa- 
rated by a distance d where channel j contains Nj par- 
ticles [33]. The total number of particles in the system 
is N = ^j^iNj. The particle-particle interactions are 
repulsive, and each particle interacts both with particles 
in the same channel as well as with particles in nearby 
channels. The number of randomly placed pinning sites 
in channel j is given by rij. All of the particles in the 
system are subjected to a uniform external driving force 
F^ . We apply periodic boundary conditions in the x di- 
rection, along the length of the channels, and the total 
length of our system is L = 72ao, where oq is our unit 
of length that is typically a micron in colloidal systems. 
The boundaries are open along the y direction transverse 
to the channels. Fig. [T[ shows a schematic of our sys- 



tern in several different configurations. In Fig.[IJa) there 
are M = 2 channels with equal numbers of particles, 
iVi = N2, and pins, ni = 712, in each channel. Figure 
[IJb) shows a two channel system with equal numbers of 
particles but with more pinning sites in the upper chan- 
nel, ni > 712, while in Fig. [He) there are equal numbers 
of pins but more particles in the upper channel, A'^i > 7V2- 
Figure [TJd) shows a multiple channel Af = 5 system with 
equal numbers of particles in each channel, Nj — Ni , and 
equal numbers of pins in each channel, rij = ni. 

The particles interact via a Yukawa or screened 
Coulomb potential that is appropriate for charged col- 
loids or for charge transport in the presence of screening. 
We show later that the same general dynamic phases also 
occur for repulsively interacting vortices in type-II super- 
conductors. For the colloidal system, the dynamics of a 
particle i is determined by integrating the following over- 
damped equation of motion: 

'^^^ Ff + F,^ + F^ + Ff . (1) 
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dt 



Here R^ is the location of particle i and 77 is the damping 
coefficient that is set equal to one. The particle-particle 
interaction force is F^'' = —^^j,i'^V{Rij)'iiij where 
V{R,j) ^ {q^E^/R,j)e^^{-nR,j), R,, = |R, - Rj|, 
Rij = (Ri - Rj)/i?y, Eq = Z*^/47reeo, e is the sol- 
vent dielectric constant, Z* is the effective charge of a 
colloidal particle, and q^ is the dimensionless squared 
colloid charge that is taken to be q^ = 2.0 unless oth- 
erwise noted. We take the screening length 1/k — 
4ao and choose the distance between channels in the 
range d < Sap (equivalent to d/n < 2) to ensure 
that particles in neighboring channels interact with each 
other. In the absence of pinning, the particles in an 
isolated channel would adopt a lattice spacing of Uj — 
Nj/L. The substrate pinning force is given by Ff — 
TZi F^{R\JR,)Q{R, ~ R\^)^ . Here particle i sits in 
channel j , rij is the number of pinning sites in channel j, 
Q is the Heaviside step function, R^ is the location of pin- 
ning site fc, RPf^ = |R, - R^l, and Rf^ = (R, - nD/R^^. the channels, V = J^^i ^^ 
The pinning site radius is Rp and the maximum strength 
of the pinning sites in channel j is F^. Within each 
channel the pins are placed in randomly chosen non- 
overlapping positions that are different for each channel. 
This pinning model has been used in previous simula- 
tions to represent quenched disorder for colloidal sys- 
tems, classical electron systems, and vortices in super- 
conductors. The external drive F-^ = F^x is applied 
in the positive x-direction. The drive is increased in 
increments of 6F^ = 0.001, and after each increment 
F^ is held fixed for 10^ simulation time steps to en- 
sure that any transient effects have subsided. The ther- 
mal force Ff is modeled by adding a Langevin noise 
term with the following properties: {Ff{t)) — and 
(Ff{t)F[{t')} = 2r]kBS,jSit - 1'), where fcs is the Boltz- 
mann constant. We initially consider the T = case. 
We measure the average particle velocity in each chan- 
nel for every force increment, Vj — {N^ 12i=i '^■i.)^ where 



FIG. 2: The velocity of each channel, Vi (dark line) and 
V2 (light red line), in a two channel (M = 2) system versus 
applied drive F with fixed quenched disorder for varied dis- 
tance d between the channels. Each channel has the same 
number of particles, A^i = N2, and pinning sites ni = n2. 
Here ni/iVi = 0.25, Ff = F^ = 6.0, Rp = 0.23ao, and the 
in-channel lattice constant ai = a2 — L/N\ = 1.5ao. (a) For 
d/ai = 1.133, there is an elastic depinning transition and both 
channels depin simultaneously, (b) At d/a\ = 1.47, channel 1 
depins first while channel 2 remains pinned. Once channel 2 
depins, the channels become dynamically locked. The locking 
transition is associated with negative differential conductivity 
(NDC) in Vi, as illustrated in the inset where dVi/dF° < 0.0 
at F° = 0.36. (c) For d/ai = 1.57 near F° = 0.39 there 
is a region where both channels are depinned but Vi 7^ V2, 
indicating that the channels are sliding past one another. Dy- 
namical locking occurs for F^ > 0.415. (d) At d/ai = 1.64, 
the recoupling transition is shifted to much higher F^ . In- 
set: The sum of the velocities V = Vi -I- V2 for d/ai = 1.47 
(dark line) and d/ai = 1.64 (light green line) shows the oc- 
currence of two step depinning; there is no NDC in V for the 
d/ai — 1.47 sample. 



{dRi/dt) -x. We also measure the net velocity of all 
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III. DECOUPLING TRANSITIONS FOR TWO 
LAYER SYSTEMS 

We first consider a two layer system in which we vary 
the interlayer distance d. We can effectively decrease 
the coupling between the layers by increasing d since the 
interaction between Yukawa particles becomes weaker for 
increasing interparticle distance. In Fig. [2] we plot the 
velocity force curves for channels 1 and 2 in a system with 
equal numbers of particles in each channel (TVi = A^2), in- 
channel lattice constant ai — a2 = 1.5ao, equal numbers 
of pins in each channel {ni — 712), ni/Ni — 0.25, Ff — 
F2 = 6.0, and Rp = 0.23 for varied d measured in terms 
of d/fli. For small d or strong coupling, the two channels 
depin simultaneously without any slipping between the 
channels, as shown in Fig.[2)Ja) for d/oi = 1.333. All the 



particles keep the same neighbors, indicating that the 
depinning is elastic. The value of F^ at which depinning 
first occurs is termed the critical depinning threshold, Fc- 

As d increases there is a transition to a state where 
the channels depin individually, as shown in Fig. HJb) 
for d/ai = 1.47. Here channel 1 depins at F^ = 0.32 
while channel 2 remains pinned until F^ = 0.36. Since 
both channels have equal numbers of pins, equal num- 
bers of particles, and experience the same driving force, 
the difference in depinning threshold arises due to the 
random placement of the pinning sites in the two chan- 
nels, which differs from channel to channel. As soon as 
channel 2 depins, the system transitions directly into a 
state where the moving channels are dynamically locked 
with each other and move at the same velocity, Vi — V2, 
without any slipping. As the inset in Fig. [2jb) shows, 
the dynamic locking transition coincides with a drop in 
Vi indicative of negative differential conductivity (NDC), 
where a system of particles shows a decrease in the ve- 
locity under increasing drive, dVi/dF^ < 0.0. The NDC 
occurs only in channel 1 at F^ — 0.36 when V2 jumps 
from zero to a finite value, which also corresponds to a 
peak in dVi/dF'^. 

Fig. \^c) shows the velocity force curves for d/ai ~ 
1.57, where the depinning threshold for channel 1 is 
again lower than that for channel 2. In this case, af- 
ter channel 2 depins the system does not immediately 
enter the moving locked phase. Instead, we find a re- 
gion 0.39 < F^ < 0.415 with Fi 7^ 1^2, 14 > 0, and 
V2 > 0. Here both channels are moving but at different 
average velocities, indicating that the particles in channel 
1 are sliding past the particles in channel 2. The chan- 
nels become dynamically locked for F^ > 0.415, where 
Vi = Vi- As djax is further increased, the dynamical 
locking transition shifts to larger values of F^ . For ex- 
ample, in Fig.[21[d) for djai = 1.64, the channels are only 
dynamically coupled for Fq > 0.742. 

In the inset of Fig. [Hd) we plot the overall velocity 
y = Vi -I- V2 versus F^ for the samples with d/ai = 1.47 
and d/ai = 1.64 from Fig. [2ljb) and Fig. [UJd), respec- 
tively. These are the curves that would be obtained in an 
experimental measurement of net transport rather than 
individual channel transport. For d/ai — 1.47, where 
NDC of Vi appeared in Fig. HJb), the inset of Fig. [IJd) 
shows that there is no NDC in the overall velocity V . 
This is because the decrease in Vi at F^ = 0.36 is ex- 
actly compensated by the increase in V2 , so there is no net 
drop in V at the second depinning transition. For both 
d/ai — 1.47 and d/ai — 1.64, V shows a characteristic 
step feature that indicates the presence of a two-step de- 
pinning transition. On the lower step in V, only channel 
1 is depinned, while on the upper step, both channels are 
flowing. For this particular set of parameters, we find 
only a weak signature in V of the relocking transition 
that occurs at higher drives; however, for other param- 
eters, we will show below that the onset of dynamical 
locking produces much more pronounced effects in V. 

By conducting a series of simulations we map the dy- 
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FIG. 3: Dynamic phase diagrams F^ vs d/ai for the AI — 2 
system from Fig. [2] with A''i = A'^2, ni — n2, ni/Ni — 0.25, 



Rp — 0.23ao, and ai 



a2 



1.5ao. P: pinned phase with Vi 



V2 = 0; C: coexistence phase where one channel is moving 
but the other is pinned; S: sliding phase where both channels 
are moving but Vi 7^ V2; L: dynamically locked phase where 
Vi = ^2 > 0. (a) The sample from Fig. i with F^ = f| = 6.0. 
For d/ai < 1.43 the system depins directly into the locked 
phase, while for d/ai > 1.43, dynamically induced locking 
occurs at a value of F^ that increases with increasing d/ai. 
Inset: A blowup of the main panel illustrating that the P 
phase reaches its maximum value of F^ at the transition from 
elastic to plastic depinning. (b) A sample with Ff — F2 — 3.0 
has the same features but the onsets of the S and C phases 
fall at higher values of d/ai > 1.6. 



namic phase diagram for F^ vs d/ai for the sample from 
Fig. m with F^ = Ff = 6.0 in Fig. [^a) and for a sample 
with Ff = F| = 3.0 in Fig. Mh). In the pinned (P) 
phase, both channels are pinned, V\ = V2 = 0. In the 
coexistence (C) phase, one channel is pinned while the 
other is moving. The sliding (S) phase has both channels 
moving but not locked, 14 > and V2 > but 14 7^ V2, 
while in the locked (L) phase both channels move to- 
gether without slipping, F^ = V2 > 0. In Fig. [3]^a), for 
d/oi > 1.43 the layers decouple at depinning and the 
width of region S grows with increasing djax. The width 
of the coexistence phase saturates as djax increases due 
to the decoupling of the particles, which causes the de- 
pinning thresholds of each channel that mark the borders 
of the C phase to be determined only by the quenched 
disorder configuration and the interactions among parti- 
cles within the channel, and to be insensitive to the po- 
sitions of the particles in the neighboring channel. The 
recoupling that marks the transition between the S and 
L phases rapidly shifts to higher F^ with increasing dj a\ 



for djax > 1.613. In Fig. ^h) we show that in a sample 
with weaker pinning of Ff = F2 = 3.0, the transition 
directly from P to L persists out to the higher value of 
d/ax = 1.6 but that the same general features found for 
higher F^ in Fig.[3I^a) still occur and are shifted to higher 
d/ax and lower F^ . 

The dynamical phase diagrams in Fig. |3] have many 
similarities to the dynamical phase diagrams obtained for 



2D and 3D driven vortex systems [IJ, Ua M |26|, |37| . The 
2D vortex system exhibits pinned, plastic, and dynam- 
ically ordered phases, with the vortex lattice ordering 
into either a moving smectic or moving crystal configu- 
ration [13--15;j . Our system has the same pinned phase, 
while the coexistence phase that contains a mixture of 
moving and pinned channels corresponds to the plastic 
flow regime for the vortex system. The sliding phase we 
observe resembles the moving smectic phase found for 
vortices, while our locked phase is equivalent to the 2D 
dynamically reordered phase or, for 3D vortex systems, 
to the transition from decoupled layers to coupled 3D 
vortex lines [23, |2^, |40] ■ In the vortex system, as the pin- 
ning strength increases or the vortex- vortex interactions 
are weakened, the reordering transition shifts to higher 
drives [14| . The dynamical reordering occurs because 
the pinning is effectively weakened when the vortices are 
moving rapidly. In our system, in the locked phase the 
effectiveness of the quenched disorder is dramatically re- 
duced. As the interchannel particle-particle interaction 
strength decreases for increasing d/ai, the quenched dis- 
order becomes more effective, the recoupling transition 
shifts to higher F^, and the additional C and S phases 
appear between the P and L phases. 

In the inset of Fig. ^a) we plot a blowup of the re- 
gion near d/ai = 1.4 for the Ff = Ff = 6.0 sample 
where the transition between elastic depinning directly 
into the locked phase and plastic depinning into the co- 
existence phase occurs. The critical depinning threshold 
Fc reaches a peak at d/ai = 1.43 at the transition, while 
for d/ai > 1.43, Fc decreases and then saturates with 
increasing d/ai. The Fc peak has all the hallmark fea- 
tures of the peak effect phenomenon found for vortices 
and other systems where a peak in the critical depinning 
threshold occurs near the transition from elastic to plas- 
tic depinning [3, [III, UM, [la|- In many of these studies 
the Fc peak occurs inside the plastic depinning regime 
just beyond the transition point, while Fc for the elastic 
depinning is always lower than for the plastic depinning 
[5| . In our two layer system, the Fc peak falls right at the 
transition out of the elastic depinning regime, and Fc de- 
creases with increasing d/oi within the plastic depinning 
regime. 

Determining whether the effectiveness of the pinning 
in the plastic or decoupled phases is higher than in the 
elastic phase can be ambiguous when different quanti- 
ties are measured [^, llM, ■ If we consider only the overall 
velocity ^ at a fixed value of F^ as we vary d/ai, we 
find that V drops in the decoupled regime, suggesting 
that the pinning is more effective in the decoupled state. 
This contrasts with the behavior of the depinning thresh- 
old Fc, which drops in the decoupled regime, suggesting 
that the pinning is less effective in the decoupled state. 
In Fig. [U we plot V versus d/ai for the sample from 
Fig. ^a) with Ff = F| = 6.0 for fixed F^ values rang- 
ing from F° = 0.4 to F° = 1.0. For F° < 0.8, V 
drops or becomes zero at the decoupling transition, and 
for 0.6 < F^ < 0.8 the value of V for d/ai above the de- 
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FIG. 4: The summed velocity V = Vi + V2 for the M = 2 
system in Fig. [3ja) with Ni — N2, ni = n2, n\/N\ = 0.25, 



Rp — 0.23ao, ai = 1.5ao, and _Ff 



Fi 



6.0 at F 
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0.4, 



0.5, 0.6, 0.7, 0.8, 0.9, and 1.0, from bottom to top. The value 
of V in the decoupled regime is always lower than in the 
coupled regime, and the dip in V at the decoupling transition 
gradually fades away as F^ increases. 
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FIG. 5: (a) The critical depinning threshold Fc vs d/ai for 
the M = 2 system in Fig. \T[a.) with iVi = A'"2, r^i = n2, 
ni/Ni = 0.25, Rp = 0.23ao, and ai — 1.5ao for varied Ff = 
Ff = 15.0, 12.0, 9.0, 6.0, and 3.0, from top to bottom. In 
each case the peak in Fc falls at the transition from elastic to 
plastic depinning. (b) Fp^, the value of Ff and F2 at which 
a transition from elastic to plastic depinning occurs, vs d/ai 
in the same system. 



coupling transition is smaller than the value of V for d/ai 
below the decoupling transition. For F^ > 0.8, the size 
of the dip in V at the decoupling transition gradually di- 
minishes. This behavior is very similar to that observed 
across the peak effect for vortices in type-II supercon- 
ductors, and indicates that for the moving particles, the 
pinning is more effective in the decoupled phase than in 
the coupled phase. The behavior of V differs from that 
of Fc because V is associated with moving particles while 
Fc is determined by the static configuration of the pinned 
particles. 

In Fig. [5|Ja) we plot Fc versus d/ai for a range of Fp 
values in the system from Fig. ^a) to illustrate that the 




FIG. 6: Vi (dark line) and V2 (light red line) vs F° for 
the M — 2 system in Fig. [SJa) with A'^i = A''2, ni = 712, 



ni/A^i = 0.25, Rp = 0.23ao, ai 



1.5ao, and Ff = Ff 



15.0 

for (a) d/ai = 1.147, (b) d/ai = 1.2, (c) d/ai = L27, and 
(d) d/ai = 1.47. Inset of (d): The dVi/dF° and dV2/dF° vs 
_F^ curves both peak at the dynamical coupling transition at 
F^ = 2.0; for _F^ > 2.0 above the transition, the fluctuations 
in both curves show increased correlation. 



decoupling transition is associated with a peak in Fc. 
On average, Fc increases with increasing F'p and with in- 
creasing d/ ai. The decoupling transition, and with it the 
peak in F^, shifts to lower values of d/ai with increasing 
F'P and becomes more prominent. In all cases the peak 
in Ff. falls exactly at the transition from elastic to plastic 
flow. For a given value of d/ai, it is possible to induce a 
decoupling transition by increasing F^] the value of Ff 
and F2 at the decoupling transition is termed F^^ . In 
Fig. [Sfb) we plot F^^ versus d/oi, where we find that 
as d/ ai decreases and the coupling between the channels 
increases, F^^ increases faster than linearly. 

The velocity force signatures illustrated in Fig. [5] are 
generally robust against changes in pinning strength. For 
samples with stronger pinning, Ff = F2 > 6.0, many 
of the features become more prominent. For sufficiently 
strong FP, the system can depin directly into the slid- 
ing state with no intermediate coexistence state, as illus- 
trated in|ni;a) for Ff = if = 15.0 and d/ai = 1.147. 
Here both channels depin simultaneously at F^ — 0.915, 
and for 0.915 < F^ < 1.025 the channels move at dif- 
ferent velocities with Vi ^ V2. We also find a larger 
number of jumps in the velocity force curves for higher 
FP, as shown in Fig. [6ja) near F^ = 1.08. The jumps 
are indicative of the solitonlike nature of the particle mo- 
tion in channels with strong pinning below the recoupling 
transition. Due to the random spatial placement of the 
pinning sites, the individual pins vary in their effective- 
ness at trapping particles. As the driving force is in- 
creased below depinning, particles escape from the least 
effective pins only to pile up behind particles that are 



trapped in the more effective pins. Since the channel is 
ID, flowing particles cannot pass trapped particles. As 
the driving force further increases, particles in the pileup 
regions approach each other more closely and exert a 
greater force on the trapped particle that produced the 
pileup. At depinning, the trapped particle escapes from 
the pin and is immediately replaced by another particle 
that becomes pinned in its place. The trapped parti- 
cle travels across the system and rejoins the pileup; the 
extra force it exerts on the particles ahead of it causes 
the trapped particle at the front of the pileup to depin 
and repeat the process. This picture is oversimplified; in 
actuality, each channel contains multiple trapping sites 
that create pileups and particles jump from one pileup 
to another above depinning. Each of the trapping sites is 
associated with some local depinning threshold F^f"^. As 
F^ increases, one trapping site after another reaches the 
condition F^""^ < F^ and ceases to trap particles. The 
resulting enhancement of the mobility of all the particles 
in the channel manifests itself as a jump in channel ve- 
locity Vj. As FP increases, there is a greater spread in 
Fjf^, producing a larger number of steps in Vi and V2 
above depinning. We note that such velocity force jumps 
are often called switching events in sliding charge den- 
sity wave systems [11]. For vortex systems, near the peak 
effect regime a series of jumps and dips in the current- 
voltage curves and dV/dl curves can appear. These have 
been termed a fingerprint phenomenon since the same 
features repeat upon cycling the IV curves, indicating 
that the features result from the detailed configuration 
of the pinning sites pll . |4^ . The curves in Fig. HJa) are 
obtained for a specific pinning realization. For a different 
pinning realization, the general features of the velocity 
force curves remain the same but the location and height 
of the velocity jumps will change. In Fig. IH^b) , a sample 
with the same disorder configuration and strength but 
with d/ai = 1.2 has separate depinning thresholds for 
channels 1 and 2 along with a sharp transition into the 
locked phase at F^ ~ 1.2. In the d/ai — 1.27 sample of 
Fig. [6l^c), we observe a strong negative differential con- 
ductivity signature in Vi at the transition between the 
sliding and locked phases. There is also a change in the 
rate at which Vi and V2 increase with F^ in the locked 
phase, with a more rapid increase for 1.25 < F^ < 1.365 
and a slower increase for F^ > 1.365. For d/a = 1.47 in 
Fig. E^d), the channels do not lock until F^ = 2.03, and 
Vi shows additional structure within the sliding phase. 
For increasing FP beyond what we show here, we find 
similar features in the velocity force curves; however, the 
number of jumps and anomalies in the curves further in- 
creases as the microscopic configurations of the pinning 
sites begin to dominate the behavior completely. 

We next compare the velocity fluctuations in the indi- 
vidual channels in the unlocked and locked state. Fig. [7| 
shows the instantaneous velocities Vi{t) and V2{t) for 
channels 1 and 2 as a function of time measured in sim- 
ulation steps for the system in Fig. ^c) at d/ai = 1.57 
6.0. At F^ = 1.395 in Fig. Ela), the 
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FIG. 7: (a) The instantaneous velocity of each channel Vi(t) 
and V2{t) versus time in simulation steps for channels 1 (dark 
black line) and 2 (light red line) for the M = 2 system in 
FigMc) with iVi = iV2, ni = 712, ni/Ni = 0.25, Rp = 0.23ao, 
ai = 1.5ao, d/ai = 1.57, and F^ = F| = 6.0 at F^ = 0.395. 
The system is in the S phase and both channels are moving at 
different velocities, (b) S'(tj) for the time series in panel (a) 
shows that channel 2 has greater spectral power at the lowest 
frequencies compared to channel 1. lj is reported in inverse 
simulation time steps, (c) Vi{t) and V2(f) for the same system 
with F^ = 0.47 where the channels are locked, (d) S{lo) for 
the time series in panel (c) shows nearly equal spectral weight 
in each channel at the fundamental frequencies. 
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FIG. 8: The dynamic phase diagram F 
M = 2 system in Fig. [2] with A^i = N2, m 
0.25, R-p = 0.23ao, ai = 1.5ao, Fl = and Fl = 15.0. P: 
pinned phase; C: coexistence phase where only one channel is 
moving; S: sliding phase where Vi 7^ V2 > 0; L: locked phase. 
As d/ai increases, the depinning threshold for channel 1 goes 
to zero and the system immediately enters region C for small 
nonzero F^ . Inset: A blowup of the main panel near the 
transition from elastic to plastic depinning. 



IV. VARIED RELATIVE DISORDER 
STRENGTH AND HYSTERESIS 



system is in the S phase and both channels are moving 
with different average velocities Vl > V2 > 0. Both chan- 
nels show a periodic response but channel 1 has a higher 
frequency and higher average velocity, while V2(t) drops 
nearly to zero during each cycle. In Fig. Eljb) we plot the 
corresponding power spectra S(ui) of the velocity signals, 
with w measured in inverse simulation time steps. Chan- 
nel 2 has greater spectral weight at the low frequency 
oi uj — 1.1 X 10"'* corresponding to its fundamental fre- 
quency, while the peak at this frequency for channel 1 is 
orders of magnitude smaller in height. Channel 1 shows a 
response at this frequency due to its coupling with chan- 
nel 2. Near uj = 1.62 x 10^**, the fundamental frequency 
of the faster moving channel 1, the peak in S{uj) is higher 
for channel 1 than for channel 2. This result shows that 
the two channels each have a periodic velocity signal cor- 
responding to a different washboard frequency; however, 
since the channels are interacting, both velocity signals 
contain both washboard frequencies. In Fig. [Tfc) we plot 
Vi{t) and V2{t) for the same system at F^ = 0.47 in 
the locked regime where the channels have the same av- 
erage velocity. Here both channels exhibit the same fun- 
damental frequency even though the shapes of Vi (t) and 
V2{t) do not completely overlap. The corresponding S{uj) 
in Fig. [Zl^d) shows that both channels have nearly the 
same spectral weight at the fundamental frequency of 
UJ — 1.48 X 10~* and its harmonics, unlike the unequal 
peaks which appeared in the unlocked phase in Fig.[71[b). 



We next consider the effects of varying the pinning 
strength from channel to channel in order to create one 
channel with high pinning strength and one channel with 
low pinning strength. Figure [SJa) shows the dynamic 
phase diagram of F^ versus d/ai for a two channel sys- 
tem with Ff = 15.0 and F2 = 0, while the inset il- 
lustrates a blowup of the area near the onset of plastic 
depinning. The strongly pinned particles in channel 1 
are able to effectively pin the particles in channel 2 via 
particle-particle interactions alone, allowing the pinned 
regime to persist for finite values of F^ in spite of the fact 
that the pins in channel 2 have zero strength. For low 
d/ai when the interactions between particles in neigh- 
boring channels are strong, the depinning of both chan- 
nels occurs simultaneously in a single step at Fc from 
the pinned state to the moving locked state. As d/ai 
increases, Fc increases and reaches a maximum value 
at d/ai — 1.16 corresponding to the decoupling tran- 
sition. For d/ai above the decoupling transition, the de- 
pinning force Fc for channel 2 decreases monotonically 
with increasing d/ai, unlike the saturation that occurred 
in Fig. [3] for samples with i^f = Fj*. As d/a 1 increases, 
the coupling between the channels decreases monoton- 
ically, and since the particles in channel 2 are pinned 
only due to their interactions with the pinned particles 
in channel 1, Fc for channel 2 decreases with increasing 
d/ai. The transition line between the C and S phases, 
which corresponds to the depinning force Fc for channel 
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FIG. 9: Vi (dark line) and V2 (light red line) vs F° under a 
cycled drive for the M = 2 system in Fig. [8] with TVi — N2, 
m = n-i, ni/Ni = 0.25, Rp = 0.23ao, oi = 1.5ao, Ff = 15.0 
and F2 — 0.0. (a) At d/ai — 1.27 we find hysteresis in 
Vi near the transition to the pinned phase but there is no 
hysteresis in V2. The arrows indicate the V2 curves obtained 
for sweeping F^ up and down, (b) For the weaker coupling 
of d/ai = 1.33, there is no hysteresis and dynamic coupling 
occurs at F° = 1.43. 



1, grows more rapidly above the decoupling transition 
with increasing d/ai in Fig. [8] than it did for samples 
with Ff = F2 in Fig. |31 This is because the particles 
in channel 2 produce relatively little effective drag for 
the particles in channel 1 when there is no pinning in 
channel 2, and this drag rapidly becomes nearly zero as 
d/ai increases and the coupling between the channels de- 
creases. The transition from the S to L phases also grows 
more rapidly with increasing d/ai in Fig. [5] than for the 
samples with Ff = F2. Our results show that coupling- 
decoupling transitions can persist even when the pinning 
in one channel is completely absent. 

For the parameters examined so far in this work we 
observe only weak hysteresis in the velocity-force curves 
occurring near the sharp jumps in Vi or V2. Hysteresis 
is often absent in ID systems, as shown with the no- 
passing rule [ijl; however, in studies using phase-field 
modeling, hysteresis has been observed in ID systems 
[6|. In our system we never find hysteresis when the 
depinning is elastic, but hysteresis can occur for plastic 
depinning, as illustrated in Fig. EJa), for a system with 
the same parameters as in Fig. |8]at d/ai — 1.27. Here 
hysteresis occurs at the depinning transition only for the 
pin- free channel 2 but not for channel 1. The summed 
velocity V — V1+V2 also exhibits hysteresis at depinning, 
but by examining the separate velocity signals we can 
determine that the hysteresis originates from only one 
of the channels. Figure [Hl^a) also shows that there is no 
hysteresis across the S-L transition at F^ = 0.526 in 
spite of the discontinuous jump in both Vi and V2 that 
occurs at this transition. In Fig. ^h) we plot Vi and 
V2 for cycled F^ in a more weakly coupled sample with 
d/ai = 1.33, further from the value of d/ai = 1.16 where 
the transition from elastic to plastic depinning occurs. In 
this case we find dynamic locking of the channels above 
F^ — 1.43; however, there is no hysteresis across any of 
the transitions. 

We can also examine the effects of varying the relative 




FIG. 10; Dynamic phase diagram F° vs Fl/F.f for an M = 2 
system with Ni — N2, ni = n2, ni/Ni = 0.182, Rp = 0.35ao, 
ai = 1.5ao, d/ai = 1.7, and fixed Fj' — 0.75. P: pinned 
phase; C: coexistence phase; S: sliding phase; L: locked phase. 
Near Ff /F2 = 1.0, where the pinning strength is the same in 
both channels, the system depins elastically. The depinning 
becomes plastic once Ff /Ff moves sufficiently far away from 
1 in either direction. 



pinning strength in the two channels by altering Ff /Ff 
for a system with iVi = N2, ni = n2, ni/Ni = 0.182, 
Rp = 0.35ao, and d/ai = 1.7 as shown in Fig. [101 Here 
F2 = 0.75 and Ff is varied. For this particular set of 



parameters at Ff /Ff = 1.0 the depinning is elastic. For 
< Ff/Ff < 0.6, the system depins plastically into the 
coexistence state where channel 1 is moving and chan- 
nel 2 is pinned, and for higher F^ the system enters the 
locking regime without passing through the sliding phase. 
At Ff/Ff = the depinning threshold of channel 1 is 
nonzero due to the strong interaction of the particles in 
channel 1 with the pinned particles in channel 2. As Ff 
increases from zero the width of the pinned region grows 
as the depinning threshold of channel 1 rises; however, for 
< Ff/Ff < 0.6 the depinning threshold of channel 2 re- 
mains constant since Ff is fixed. For 0.6 < Ff/Ff < 1.0 
the depinning is elastic and both channels depin simul- 
taneously; further, within this regime the pinned region 
has the most rapid growth for increasing Ff/Ff. For 
Ff/Ff > 1.0 the depinning becomes plastic again and 
each channel depins separately; however, channel 2 now 
depins first since Ff > Ff. As Ff/Ff increases above 1, 
the width of the pinned phase passes through a small dip 
but remains nearly constant since the depinning thresh- 
old of channel 2 is determined by Ff which is held fixed. 
Near Ff/Ff = 1.2 we find the onset of the sliding phase, 
which grows rapidly in width with increasing Ff/Ff. 
The sliding phase is followed at higher F^ by the dy- 
namically locked phase. If we increase the fixed value 
of Ff and sweep Ff/Ff, the extent of the elastic depin- 
ning window diminishes as Ff becomes larger until for 
sufficiently large Ff there is no longer an elastic depin- 




FIG. 11: Dynamic phase diagram F^ vs Ff for an M = 2 
sample with F| = Ff, Ni ^ N2, ni = n2, ni/Ni = 0.25, 
Rp — 0.23ao, ai — 1.5ao, and d/ai = 1.7. P: pinned phase; 
C; coexistence phase; S: shding phase; L: locked phase. The 
depinning is elastic for Fp < 0.8. Inset: dFc/dFf for the 
same system shows a pronounced dip near the transition from 
elastic to plastic depinning. 



FIG. 12: Dynamic phase diagram _F^ vs density 1/ai for an 
M — 2 sample with A''i — N2, ni — n2, ai =02, ni/Ni — 
0.25, Rp = 0.35ao, d = 2.75ao, and Ff = F| = 3.0. P: pinned 
phase; C: coexistence phase; S: sliding phase; L: locked phase. 
For 0.07 < 1/ai < 0.305, the system transitions from a pinned 
phase into a pulsed coexistence phase PC. 



ning regime. Instead, the sliding phase appears for all 
values of Ff /F2. We have also considered samples with 
F^ = F2 but with different pinning densities ni 7^ n2, 
and we find the same general features as a function of 
ni/n2 as we have shown here for varying Ff / F2. 

In Fig. [TT] we plot the dynamic phase diagram for F^ 
versus Ff in a sample with equal pinning strength in both 
channels, F| = Ff , and with d/ai = 1.7. For Ff < 0.7, 
the depinning is elastic and the depinning force Fc drops 
to zero as Ff drops to zero. For Ff > 0.8, the depinning 
is plastic and all of the dynamical transition lines shift to 
higher F^ for increasing Ff , with the S to L transition 
increasing the most rapidly. In the inset of Fig. [TT] we 
plot dFc/dFf versus Ff where we find a roughly linear 
increase in the regime < Ff < 0.55. This implies that 
Fc oc (Ff)^, which is consistent with elastic depinning 
in the collective regime [IJ]. Just above the transition 
to plastic depinning, dFc/dFf passes through a strong 
drop and then recovers to a nearly constant value near 1, 
consistent with Fc ex F^ as expected for the single parti- 
cle limit. This results indicates that the transition from 
elastic to plastic depinning is associated with clearly ob- 
servable features in the transport curves. If we fix Ff and 
instead increase the density of pinning sites while hold- 
ing ni — n2, we find dynamical phases with tendencies 
similar to those shown in Fig. [TT] 



V. DENSITY DEPENDENCE AND DYNAMIC 
COMMENSURATION EFFECTS 

We next consider the effect of holding the distance d 
between channels fixed while varying the density of the 



particles which changes the lattice constants ai and 02. 
In Fig. [T2]we plot the dynamic phase diagram as a func- 
tion of F^ and density 1/ai for varied ai = L/Ni in a 
system with an equal number of particles in each channel, 
Ni = N2, Ff = Ff = 3.0, and with fixed d = 2.75. The 
depinning threshold Fc generally decreases with increas- 
ing density due to the increase in the particle-particle in- 
teractions relative to the pinning strength as the particles 
get closer together within each channel. For intermediate 
densities of 0.305 < 1/ai < 0.47, the depinning is elastic, 
while for 0.07 < l/oi < 0.305, the depinning is plastic 
and the sample enters either the coexistence phase C, the 
shding phase S, or a pulsed coexistence phase PC. Each 
channel is alternately pinned or moving in the PC phase, 
so that this phase resembles the C phase except that the 
stationary and moving channels keep switching places. 
At lower densities 1/ai < 0.305, the system no longer 
depins elastically due to the increase of the effectiveness 
of the pinning, which induces a decoupling of the chan- 
nels. For l/fli < 0.07, we have A^i > ni. Since there are 
more pinning sites than particles, the depinning occurs 
in the single particle limit and it is no longer meaning- 
ful to describe the pinning as plastic or elastic since the 
particles in the channels are essentially noninteracting at 
depinning. Additionally, since Ff and Ff are held fixed, 
at these low densities all the particles in each channel 
depin simultaneously and move at the same velocity so 
the C and S phases cannot occur. 

Figure [T^ shows that for 1/ai > 0.47, the depinning is 
plastic again and the S-L transition shifts to higher F^ 
with increasing 1/ai. This might seem surprising since 
at higher densities the interactions between particles in 
the same channel become stronger; however, it is known 
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from studies of layered vortex systems that increasing the 
field can reduce the coupling between layers for fields well 
below Hc2 [401 ■ ^^ the density increases, the interactions 
between particles in the same channel become stronger 
more rapidly than the interactions between particles in 
neighboring channels. This enhances the S phase since 
it becomes more difficult for the channels to lock. In or- 
der to minimize energy, the moving particles adopt a zig 
zag configuration with the particles in one channel offset 
relative to the particles in the other channel. This ap- 
proximates a triangular lattice configuration. As 1/ai in- 
creases, the lattice spacing along the channels shrinks but 
the lattice spacing perpendicular to the channels is fixed, 
producing an effectively increasingly anisotropic triangu- 
lar lattice configuration. Eventually the spacing becomes 
so anisotropic that there is very little energy difference 
between an alternating configuration and one in which 
the particles are adjacent to each other, so the channels 
decouple. 

The phase diagram in Fig. [T^] exhibits a number of 
features found in phase diagrams for strongly layered su- 
perconductors containing disordered pinning sites. For 
example, the superconducting system undergoes a low 
field decoupling transition when the pancake vortices in 
each layer are far apart and only weakly interacting, while 
at higher fields there is a field induced decoupling [40| . 
In experiments, samples show both a low field disorder- 
ing transition and a high field decoupling transition [4J| . 
It was argued that the low field transition is pinning in- 
duced and occurs when the pinning energy overwhelms 
the weak vortex- vortex interactions at low vortex density, 
while the higher field disordering transition indicates the 
onset of the peak effect [4J]. In our system, as we in- 
crease PP the lower and upper endpoints of the elastic 
depinning regime approach each other until the elastic 
depinning disappears completely. Conversely as we de- 
crease PP the elastic depinning region becomes more ex- 
tended. Although there is a small increase in Fc at the 
lower transition from plastic to elastic depinning, we find 
no feature in Fc at the upper transition from elastic to 
plastic depinning; however, if we measure the average 
particle velocity for fixed drive, the velocity drops across 
the higher density transition, similar to the effect illus- 
trated in Fig. m 

We next examine the effects of varied particle ratio for 
a sample with Ff — F2 — 6.0, ni — n2, d/a2 ~ 1.33, and 
712 /7V2 = 0.167. We fix N2 and vary Ni. In Fig. [TSlJa) 
we plot 2Vi/N and 2V2/N versus F^ for a sample with 
N1/N2 = 1/2. Here the velocities are normalized by N/2, 
where N = Ni + N2. In this case, channel 2 contains a 
larger number of particles and depins at lower F^ than 
channel 1. The slope of 2V2/N is also greater than that 
of 2Vi/N. We replot the same curves in Fig. [THTb) with 
the normalizations V{ — Vi/Ni and V2 — V2/N2- Under 
this normalization the curves have almost the same form 
as the velocity force curves in samples with an equal num- 
ber of particles in each channel, and it is now possible to 
distinguish the transition into a locked phase, which also 




FIG. 13: Channel velocities vs F^ for an M = 2 system with 
F[ = F| = 6.0, m =712, n2/N2 = 0.167, Rp = 0.23ao, and 
d/a2 = 1.33 for fixed N2 and varied iVi. (a) 2Vi/N (lower 
dark line) and 2V2/N (upper light red line) for N1/N2 = 
1/2. The velocities are normalized by N/2, where the total 
particle number N = Ni + N2. Channel 2 contains a larger 
number of particles than channel 1 and therefore V2 increases 
more rapidly than Vi. (b) The same curves plotted as Vi = 
Vi/Ni (lower dark line) and V2 ~ V2/N2 (upper light red 
line). Here a clear dynamical locking transition appears. Also 
shown are the corresponding dVi /dF^ (heavy green line) and 
dV2 /dF^ (light blue line) curves, (c) Vi (lower dark line) 
and V2 (upper light red line) for N1/N2 = 0.375. Here no 
dynamical locking occurs within this range of F^ . (d) V{ 
(lower dark line) and V2 (upper light red line) for N1/N2 = 
0.292 where dynamical locking occurs. 



appears as a signature in the dV^ /dF^ and dV^ /dF^ 
curves shown in Fig. [TSf b). At other fillings such as 
N1/N2 — 0.375 shown in Fig. fTST c). there is no dynam- 
ical locking within this range of F^ . For lower fillings, 
dynamical locking appears again as shown in Fig. llSf d) 
for N1/N2 — 0.292. We find dynamical locking only for 
N1/N2 = 1.0, N1/N2 = 1/2, and N1/N2 < 0.333. 

Fig. [HTa) shows the dynamic phase diagram for F^ 
versus N1/N2 for the system in Fig. [T31 The solid line 
at N1/N2 — 1.0 indicates that at this ratio the depin- 
ning is elastic and the system passes directly from the 
pinned to the moving locked phase. For N1/N2 < 1.0, 
channel 2 depins first since it contains a larger number 
of particles, and the depinning threshold of channel 1 
grows with decreasing N1/N2 until reaching a maximum 
value of Fc = Ff = 6.0 at the lowest values of Ni/N2. 
The dotted fine at N1/N2 = 1/2 in Fig. [HJa) indicates 
a transition from the S to the L phase at P^ = 4.5 as 
illustrated in Fig.fnjb). For lower fillings N1/N2 < 0.33 
the system enters the locked state above P^ = 7.0. For 
N1/N2 > 1.0, channel 1 depins first since it now contains 
a larger number of particles than channel 2, and the de- 
pinning threshold for channel 2 remains nearly constant 
as N1/N2 increases. We note that there could be a dy- 
namical locking at external drives much higher than the 



11 



Ph4 
2 
0, 



-l' ! ' 


' ' '(a)- 




- 


A ' 


- 


^X>=^ 


\ 4- 1 1 1 = 



0.4 

°B.0.3 

0.2 



1 ' 
c 


f^ _^ 


- ^ 1 


^ c : 


P 


I 


" , 1 , 


1 , i> 



0.5 1 1.5 



0.5 1 1.5 




FIG. 14: (a) Dynamic phase diagram F^ vs N1/N2 for the 
M = 2 system in Fig. [13] with m = 712, n2/A^2 = 0.167, i?p = 
0.23ao, d/a2 = 1.33 and Ff = Jf = 6.0. P: pinned phase; C: 
coexistence phase; S: shding phase; L: locked phase. The line 
at N1/N2 — 1.0 indicates that for this filling only the system 
depins elastically directly into the L phase. The dashed line at 
N1/N2 — 0.5 indicates a transition from the S phase into the 
L phase at this filling, (b) A blowup of panel (a) in the region 
near the depinning line shows dips in the depinning force Fc 
at N1/N2 — 0.5 and 1.0. (c) Dynamic phase diagram F^ 
vs q^, where q^ is the dimensionless squared colloid charge, 
for a sample with rii = n2, n2/N2 = 0.167, N1/N2 = 1, 



Rp = 0.23ao, Ff 



F! 



6.0, and d/ai 



1.33. PC: the 



pulsed coexistence phase described in Fig. 1121 For strong 
interactions (high q^) the depinning is elastic. For small q^ 
the depinning threshold is high. Inset: dFc/dq^ vs q^ for the 
same sample shows that at q^ = 1.15 at the transition from 
plastic to elastic depinning there is a peak corresponding to 
a change in the slope of the depinning curve. 



range of F^ wc consider here. 

In Fig. iMT b') we plot a blowup of the phase diagram 
from Fig. I14f a) near the depinning transition. Clear local 
minima in Fc appear at N1/N2 = 1.0 where the depin- 
ning is elastic and at N1/N2 = 0.5 where the depinning 
is not elastic but where the dynamical locking occurs at 
a much lower drive than for nearby fillings. The behavior 
at these two fillings can be viewed as a commensuration 
effect, where the coupling between the layers is enhanced 
at integer and certain fractional ratios of the filling fac- 
tors. Commensuration effects have been studied in a 
variety of solid-on-solid systems where there is match- 
ing between the number of particles and the number of 
potential minima [45]. An example of such a system is 
vortices in type-II superconductors with periodic pinning 
sites where peaks in the critical depinning force occur at 
integer [4^] and rational [43] matching fields. In our sys- 
tem commensuration occurs not between the number of 
particles and the number of pinning sites but between 
the number of particles in the different channels. The 
depinning force is reduced at the commensurate fillings 
due to the enhanced coupling between the layers. The 



strongest dip in Fc occurs when the depinning is elastic 
at N1/N2 = 1.0, while for fillings just above and be- 
low this value Fc passes through local maxima. At the 
incommensurate fillings, there are geometrically neces- 
sary topological defects that enter the zig-zag structure 
formed by the particles in the two channels. These de- 
fects effectively soften the structure and reduce the cou- 
pling between the layers. This is similar to what occurs 
in the peak effect where the depinning force is higher for 
plastic depinning than for elastic depinning. We find no 
apparent anomaly at N1/N2 = 2.0 in Fig.fHTb): however, 
we expect that for other parameter values, more commen- 
suration effects would be observable and that additional 
dynamical locking regimes would also appear at other 
fillings. Commensuration effects are also observed in two 
layer and three layer systems without pinning when only 
one layer is driven [39|. In these systems the commen- 
suration effects occur at fillings where the channels are 
more strongly coupled, so that the drive at which rela- 
tive slip begins to occur between the channels is much 
higher than for incommensurate fillings. Our results in- 
dicate that a distinct commensuration effect can occur in 
layered systems that differs from commensuration effects 
observed for particles moving over fixed substrates. 



It is also possible to change the overall particle-particle 
interaction strength by varying q^ . This can be achieved 
for colloidal systems by altering the effective charge on 
the particles or by changing the screening length. For 
vortex systems the vortex-vortex interactions can also 
change significantly due to thermal effects near Tc- In 
Fig. [TDJc) we plot the dynamic phase diagram for F^ 
versus q^ for the system in Fig. [MFa) with N1/N2 ~ 1, 
Ff = Ff = 6.0 and d/as == 1.33. For g^ > 1.23 or 
strong particle-particle repulsion the system depins elas- 
tically and as q^ increases Fc gradually decreases. Un- 
like the case for constant q^ but increasing density l/oi 
shown in Fig. [T^l we find in Fig. \W[c) that there is no 
elastic-plastic transition at high q^ since the zig-zag par- 
ticle configuration is not affected by increasing q^. When 
q^ is lowered, Fc increases and a transition from elastic to 
plastic depinning occurs at q"^ — 1.15. For g^ < 1.15, the 
system depins into the pulsed coexistence (PC) phase de- 
scribed earlier in Fig [121 The transition into the L phase 
also shifts to higher F^ as q^ decreases. The overall fea- 
tures of the phase diagram in Fig. [HTc) resemble those 
of the dynamic phase diagrams constructed for 2D vor- 
tex systems when the vortex-vortex interaction strength 
is varied |15l.[48J. In the vortex case, Fc increases as the 
vortex-vortex interactions become weaker while the dy- 
namical ordering transition shifts to higher values of the 
driving force. In the inset of Fig. [TiT c) we plot dFc/dq^ 
versus g^, showing more clearly the change in the slope 
of Fc at the transition from plastic to elastic depinning 
at F^ = 1.15. 
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15: V, vs F° for an M 



8 system with Nj = Ni, 
Tij = ni, rii/Ni — 0.17, Rp — 0.23ao, ai — 1.5ao, and 
d/ai = 1.33. (a) At FJ = Ff = 3.0, the depinning is 
elastic and aU channels depin simultaneously into the locked 
state, (b) At F^ = F[ = 5.0, plastic depinning occurs and 
groups of channels lock together with each group depinning 
at a different value of F^ . At high drives all the channels 
lock. P: pinned phase; C: coexistence phase, defined as be- 
ginning when the first channel depins and ending when the 
final channel depins; L: locked phase, (c) At F^ = F[ = 7.0, 
more channels depin individually. A window of the sliding 
(S) phase appears where all channels are moving but groups 
of channels are separately locked, (d) At F? — Ff = 11.0, the 
same four dynamical phases appear but are shifted to higher 
values of F^ . 



VI. MORE THAN TWO CHANNELS 



We next consider the case of M > 2 coupled driven 
channels. Figure [15] shows the velocity-force curves Vi 
through Vs versus F^ for an M = 8 channel system with 
Nj = Ni, Uj = ni, d/ai = 1.33, and ni/Ni = 0.17. Note 
that the boundary transverse to the channels is open and 
that there is therefore no significant direct interaction 
between the particles in channel 1 and those in chan- 
nel 8. In Fig. Ulla) at Ff = Ff = 3.0, the depinning 
is elastic and the system immediately enters the locked 
phase upon depinning. All the channels depin simulta- 
neously and the velocity curves overlap completely. In 
Fig. ITSl b) at F^ = Ff = 5.0 the depinning is plastic. 
Channels 1, 2, and 3 depin first and remain locked with 
each other, followed by the depinning of a group of four 
additional locked channels and then by the depinning of 
the final channel. At F^ = 0.182 all the channels be- 
come dynamically locked together. For F^ = Ff = 7.0 
in Fig. llSf c). the system exhibits a hierarchy of dynam- 
ical locking transitions beginning at F^ = 0.186 when 
the first three channels to depin become locked together. 
This is followed by a second dynamic locking of four dif- 
ferent channels at F^ — 0.25. Channel 8 dynamically 
locks with these four channels near F^ = 0.335, and a 



FIG. 16: Dynamic phase diagram F vs Ff for the AI — 8 
system in Fig. [15] with Nj = Ni, Uj = m, Ff = Ff , ni/Ni = 
0.17, Rp = 0.23ao, ai = 1.5ao, and d/ai — 1.33. P: pinned 
phase; C: coexistence phase; S: sliding phase; L: locked phase. 




FIG. 17: Dynamic phase diagram F^ vs d/ai for the M = 8 
system in Fig. [T6l with FJ = Ff = 3.0, Nj = A^i, Uj = m, 
ni/Ni = 0.17, Rp = 0.23ao, and fixed ai = 1.5ao. P: pinned 
phase; C: coexistence phase; S: sliding phase; L: locked phase. 
Inset: Blow up of the region near the transition from elastic 
to plastic depinning shows that the depinning force F^ peaks 
at the transition. 



final locking of all the channels with each other occurs 
at F^ = 0.37. As F^ increases the number of distinct 
dynamical locking transitions increases, as illustrated in 
Fig. [T5i;d) for Ff = Ff = 11.0, while the final transi- 
tion at which all of the channels lock together shifts to 
higher values of F^ . Figure [TSl^ d) also shows that it is 
possible for individual channels to exhibit negative differ- 
ential conductivity at the onset of channel locking; how- 
ever, there is generally no NDC in the summed velocity 

Figure[Tn] shows the dynamic phase diagram F^ versus 
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FIG. 18: Vj vs F" for an M = 8 system with Nj / TVi, 
rij = m, Ff = Ff = 2.0, Uj/N'^'i^'^ = 0.167, i?p = 0.23ao, and 
d — 1.7ao. Each channel has one of two possible densities: 



d/flj 



1.33 for TV, 



iV 



high 



or d/uj 



AU of the velocities are normalized by V* = Vj/N^'^^. (a) A 
sample with alternating density Ni = N3 ^ N5 ^ N7 ^ iyhigh 
and N2 = N4 = Ne = Ng = N^°^ . All eight channels depin 
at different values of F^ and undergo a series of locking and 
unlocking transitions until reaching a state at higher drives 
where all the low density channels are locked together and all 
the high density channels are locked together, (b) A sample 
with the same parameters but with segregated density A'^i = 
N2 = N3 = N4 = iV'°" and N5 = Ne = N7 = Ns = N"^'^"^ 
shows a two step depinning process. 



Ff for the M ^8 system in Fig. [H] with Fj = Ff . The 
broad features of the phase diagram are the same as those 
found for the M = 2 case in Fig. [TT] however, there are 
multiple depinning transitions within the C phase and 
multiple locking transitions within the S phase that are 
not marked in Fig. 1161 In Fig. [T7] we plot the dynamic 
phase diagram F^ versus d/ai for the same M = 8 sys- 
tem with Fj 



Ff = 3.0 and fixed ai = LSap. At low 
d/ai the system depins elastically, and a transition to 
plastic depinning occurs near d/ai = 1.4. For d/ai > 1.4, 
the transition between the S and L phases shifts to higher 
F^ with increasing d/ai, while the general shape of the 
phase diagram resembles that of the M = 2 system shown 
in Fig. [3l The C and S phases again contain multiple de- 
pinning and dynamical locking transitions, respectively. 
We observe a larger amount of hysteresis in M > 2 sys- 
tems than in the M = 2 system. In the inset of Fig.lTTlwe 
plot a blowup of the F^ versus d/ai phase diagram near 
the depinning transition to show that Fc passes through 
a peak at the transition from elastic to plastic depinning, 
similar to the peak found for the M = 2 system in the 
inset of Fig. [3lja). 

For systems in which the number of particles per chan- 
nel is allowed to vary, Nj ^ Ni, the dynamic behavior 
becomes more complex; however, it is possible to identify 
several general features. In Fig.[TH]we plot V* versus F^ 
for an Af = 8 channel with F^ = Ff = 2.0 and d — I.Toq 
where each channel contains either A^'"sh particles with 
d/uj = 1.33 or A^'°" particles with d/a^ = 1.2267. The 
velocities are normalized by V* = Vj/N^^^^, and the 
number of pins in each channel is rij/N^^^^ = 1.67. In 
samples where all channels have equal numbers of parti- 



cles, Nj = TVi = N^'s^ or Nj = Ni = 7V'°™, the depin- 
ning is elastic for this set of parameters. The depinning 
becomes plastic when the number of particles differs from 
channel to channel. Figure [TST a) illustrates an alternat- 
ing density system with Ni = N3 = N^ = Nr = N^'<^^ 
and N2 ^ N4 ^ Ne = Ns = A^'°". Here every chan- 
nel has a unique depinning threshold. There are mul- 
tiple locking and unlocking transitions as channels be- 
come locked with a nearest neighbor channel at lower 
drives only to unlock at higher drives. The locking at 
lower drives occurs when the pinning is still able to cre- 
ate density inhomogeneities within individual channels. 
Particles form higher density regions in the areas behind 
the more effective pins where pileups occurred below the 
depinning transition. These higher density regions cre- 
ate an additional energy barrier for relative slip between 
the particles in neighboring channels, favoring the cou- 
pling of adjacent channels. This process creates a window 
0.105 < F^ < 0.115 where all of the channels couple to- 
gether, Vj — Vi. As F^ increases, the effectiveness of the 
pinning is reduced and the density within each channel 
becomes more homogeneous. At the same time, the over- 
all effectiveness of the pinning in channels with A^'"s'^ is 
reduced relative to that of the pinning in channels with 
^low j^g ^Q ^]-^g stronger particle-particle interactions in 
the higher density channels. The velocity y'^'sh jjj chan- 
nels with iV'^'sh jg therefore higher than the velocity V^°™ 
in channels with N^°^ , F'^'^h > ^low^ .^j^j^j^ ^^^^^^ jg_ 
coupling of adjacent channels in the alternating density 
system. For high enough drives F^ > 0.14, the system 
enters a doubly locked phase where all channels contain- 
ing the same number of particles are locked together, 
Vi = V3 = V5^ Vr and V2 = V4 = Ve = Vs, but Vi > V2 
by a small amount so that the two sets of locked channels 
slip with respect to each other. 

There are many other possible ways to arrange the 
Af = 8 channels such that half of the channels contain 
N^^^^ particles and half of the channels contain ]\f^°™ par- 
ticles. For example. Fig. [TST bl shows V* versus F^ in a 
density segregated system with Ni = N2 = N3 = N4 — 
7Vi°" and N5 = Ne = Nr = Ns = N^'^^. Here each 
group of channels with equal density acts as a unit and 
depins elastically into a locked moving state; however, 
the two groups of channels do not lock with each other. 
The depinning occurs in two steps with a lower depin- 
ning threshold for the set of channels containing N^°™ 
particles. This behavior resembles that of an iV = 2 
sample with N1/N2 slightly below 1, as in Fig. [Hfb). 
which depins into the coexistence phase and then enters 
the sliding phase at higher F^. In Fig. [TST b) we find 
that V{ > Vi^ just above the second depinning transition 
but V{ < V^ at high F^ , with a crossing of the velocity 
curves occurring near F^ = 0.105. 

In samples with an equal number of particles per chan- 
nel, particles in neighboring channels can organize into a 
two-dimensionally ordered state in order to minimize the 
particle-particle interaction energy. This ordered con- 
figuration can be perturbed by changing the number of 
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FIG. 19: The depinning thresholds _Ff "' (lower curve) and 
(upper curve) vs R/M for the M = S sample from 



Fig. ITST a) with rij 



d/aj 



1.216 for TV 



1.7ao, d/aj 



Ff = 2.0, n./AT'^'sh ^ Q^Qj^ 



1.133 for iV, 



Af" 



and 



There are R evenly spaced 



channels with Nj — N^^^^ and the remaining M — R channels 
have Nj — jv'°^. When the system is commensurate with 
Nj = Ni at R/M = or R/M = 1, the depinning is elastic. 
For other values of R/M the depinning is plastic, and the 
maximum depinning thresholds appear when geometrically 
necessary dislocations form in the system. 



particles in one or more channels, which causes defects 
to appear. To study this, we consider an M = 8 system 
in which all channels initially contain Nj = iVi°* par- 
ticles with d/uj = 1.133. We select R channels evenly 
spaced across the system with < R < M and in- 
crease the number of particles in the selected channels to 



Nj = iV'^'sii with d/aj = 1.216. As R varies, the system 
passes from a state in which all channels have Nj — 7Vi°^ 
at R/M = to one in which all channels have TV,- = TV'^'sh 
at R/M = 1. Fig. [m shows F^*^'"*, the value of F° at 
which the first channel depins, and F^'^*', the value of 
F^ at which the final channel depins, for different values 
of R/M . The geometrically necessary dislocations that 
form for intermediate values of R/M are aligned with the 
direction of the drive. At R/M = and 1, dislocations 
are absent and the depinning transition is elastic, but 
for other values of R/M when dislocations are present, 
the depinning is plastic and the depinning thresholds are 
shifted to higher values of F^ . This can also be seen in 
Fig. [THJ the alternating filling sample in Fig. [T5Fa) con- 
tains numerous topological defects and has much higher 
depinning thresholds than the segregated filling sample 
in Fig. [TST b). where there are only dislocations along the 
interface between the two fillings. This result again indi- 
cates that in our model the higher depinning forces are 
associated with plastic depinning where the channels can 
slide past each other due to the presence of topological 
defects created either by geometrical necessity for chan- 
nels with unequal particle numbers or by the quenched 



FIG. 20: The time trace of the velocities of the individ- 
ual channels, Vj{t), for the M — S sample from Fig. I18f a) 
with Ml ^ m ^ Ms ^ Nj = iV'"sh and N2 = N4 = Ne = 



Ns 



N" 



(a) At F^ 



0.1 



the system is in the coexis- 



tence phase. Channels 3 and 4 are pinned, channel 5 flows 
intermittently, channels 1 and 2 are locked at a low average 
velocity, and channels 6, 7, and 8 are locked at a higher av- 
erage velocity, (b) At F^ = 0.09 there is considerably more 
intermittency. Here channels 1, 2, 7 and 8 are locked continu- 
ously, while channels 3, 4, 5, and 6 show intermittent locking 
with the other channels, (c) At F^ — 0.125 there is less inter- 
mittency and a larger number of the channels are locked with 
other channels, (d) At F^ — 0.16, channels 1, 3, 5, and 7 
are locked and due to their higher density move with a higher 
velocity than the remaining channels 2, 4, 6, and 8, which are 
locked together. 



disorder. The results in Fig.lTSland Fig. [Tulare similar to 
previous work on mixtures of particles in the presence of 
a weak disordered substrate, which showed that the de- 
pinning force is minimized when only one particle species 
is present and there are no topological defects, but that 
for intermediate mixtures of particle species, topological 
defects appear and increase the depinning force [49|. 

The locking-unlocking transition that occurs for the 
alternating density sample in Fig. [TBJa) can be better 
characterized by analyzing the time series Vj{t) of the 
velocities of the individual channels at different values of 
F^. In Fig.Iinia), where we plot Vj{t) at F^ = 0.08, the 
system is in the coexistence regime and channels 3 and 
4 are pinned. Channel 5 switches between a temporarily 
pinned state and a sliding state, with occasional jumps 
to higher velocity flow. Channels 1 and 2 are locked to- 
gether and have an average velocity of Vi = V2 = 0.052, 
while channels 6, 7, and 8 are locked together at a slightly 
higher average velocity of Ve = V7 = Vg = 0.055 and 
show a pronounced modulation at a frequency which 
matches the frequency at which V5 drops to zero. The 
point at which V5 reaches zero corresponds to the point 
at which Vq, V7, and Vg reach their maximum values. 
Numerous periodic oscillations appear in the velocities. 
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such as the washboard frequency of the particles mov- 
ing over the pinning sites. Additionally, each channel 
has a periodicity determined by the value of Nj, and 
since this frequency differs from channel to channel, the 
channels sliding past one another experience a dynamic 
periodic potential caused by the periodicity of the parti- 
cles. This slipping process is not completely periodic but 
shows some changes over time. For lower drives, more 
channels are pinned and the velocity signals become in- 
creasingly periodic. We note that for the M = 2 case 
illustrated in Fig. [7] the velocity time traces were domi- 
nated by the fundamental frequencies which were wash- 
board signals. 

Figure I^CTb) shows Vj{t) in the same sample for F^ — 
0.09. Here there is considerable intermittency and sev- 
eral channels lock and unlock. Channels 1, 2, 7, and 8 are 
locked and move at Vi = V2 ^ Vj = Vg = 0.0625. Chan- 
nel 3 has the lowest average velocity and intermittently 
locks with channel 4, which has the second lowest average 
velocity. Channel 4 locks intermittently with both chan- 
nels 3 and 6, as well as much more infrequently locking 
with the group of channels 1, 2, 7, and 8. Channel 6 also 
locks intermittently with the channel group 1, 2, 7, and 
8 as well as with channel 4. Channel 5 has the highest 
average velocity and intermittently locks to the lower ve- 
locity value of the group of channels 1, 2, 7, and 8. In 
Fig. [20l c) at F^ — 0.125, the amount of intermittency 
is reduced. Channels 2, 4, 6, and 8 are locked together 
and have a lower average velocity, while channels 5 and 7 
lock together and have a higher but strongly fluctuating 
velocity. Channel 1, which moves independently, reaches 
the highest velocity values, while channel 3, which also 
moves independently, reaches the lowest velocity values. 
At F^ = 0.16, shown in Fig. ^d), the high density 
channels 1, 3, 5, and 7 are locked together and move at a 
higher average velocity while the low density channels 2, 
4, 6, and 8 also lock with each other and move at a lower 
average velocity. 

The multiple channel system exhibits many of the fea- 
tures predicted in mean field studies of phase slip sys- 
tems 13 7| , such as the coexistence of sliding and pinned 
phases. There are, however, many features of the mul- 
tiple channel system that are not captured by the mean 
field results, such as the multiple coupling and decoupling 
transitions between pairs of channels or groups of chan- 
nels or the peak effect in the depinning threshold at the 
transition from elastic to plastic depinning. Additionally, 
for both the two channel and multiple channel systems, 
the case where different amounts of disorder exist in each 
channel has not been studied with mean field models. We 
note that that there are many more parameters that can 
be explored for the multiple channel systems, such as 
having different densities or pinning strengths for each 
channel; however, the results for the dynamic phase dia- 
gram for the multiple channel system suggest that many 
of the generic features observed for the two level system 
will persist in the multiple channel system. We will ex- 
plore this in more detail elsewhere. 



VII. SUMMARY 

We have examined systems of two or multiple one- 
dimensional channels of coupled particles that interact 
repulsively within each channel and between the chan- 
nels. The particles are uniformly driven with an exter- 
nal drive and in the presence of quenched disorder show 
a series of dynamic phases including a pinned phase, a 
coexistence between pinned and sliding phases, a slid- 
ing phase where the channels move at different average 
velocities and slip past one another, and a dynamically 
locked phase where the particle positions in the chan- 
nels become locked and the channels move at the same 
velocity. The transition between these different phases 
can be observed as clear features in the velocity force 
curve characteristics. These features include a sudden 
drop in the velocity or the onset of negative differential 
conductivity in one of more of the channels at the dy- 
namically induced locking transition. For weak pinning 
or strong channel coupling, the depinning occurs elasti- 
cally without any slipping between channels. When the 
system parameters are varied, such as by decreasing the 
coupling between channels, we find a transition to plas- 
tic depinning, with the initial flow above depinning in the 
coexistence regime where the channels depin separately. 
At the transition between elastic and plastic depinning, 
a peak in the depinning force appears which resembles 
the peak effect phenomenon found in more complicated 
models with transitions between elastic and plastic de- 
pinning. We have also examined channels containing un- 
equal numbers of pinning sites or particles and find that 
even for a channel containing no pinning, there can be 
a finite depinning threshold due to the coupling of the 
particles in the pin-free channel with the particles in a 
channel containing pinning. For unequal numbers of par- 
ticles in the channels we observe commensurability effects 
where the depinning threshold drops when the depinning 
becomes elastic for particle number ratios at which the 
coupling between particles in neighboring channels is en- 
hanced, such as at 1:1 or 2:1. At the incommensurate fill- 
ings the channels always depin plastically. For multiple 
channels we find a hierarchy of dynamically induced cou- 
pling phases in which different groups of channels couple, 
while at higher drive all of the channels become coupled. 
Despite the apparent simplicity of our model, we find 
that even the two channel system exhibits many of the 
prominent features of depinning in more complex mod- 
els such as the existence of multiple dynamical phases as 
well as elastic to plastic depinning transitions of the type 
that have been studied for more complicated systems in- 
cluding 3D layered systems. 

The specific model studied in this paper could be re- 
alized using colloidal particles in coupled ID channels or 
in superconductors with a one-dimensional corrugation 
combined with random pinning where the vortices move 
along the easy flow direction of the corrugation. Other 
possible realizations include coupled channels of ID wires 
where Wigner crystallization may occur. 
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